General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



DE PESQUISAS ESPACIAI3 





INDEX 


Page 

ABSTRACT vi 

1 - INTRODUCTION AND APPROACH 1 

2 - INTERPRETATION METHODOLOGY AND RESULTS 3 

2.1 - VISUAL INTERPRETATION 3 

2.2 - AUTOMATIC INTERPRETATION 5 

2.3 - AREAL EVALUATION OF BARE SOIL, ARBOREAL VEGETATION AND 

PASTURE/SUGAR CANE USING LANDSAT AUTOMATIC INTERPRETATION... 10 

2.4 - CROP IDENTIFICATIONS AND AREA ESTIMATES - 1975 11 

2.5 - CROP IDENTIFICATIONS AND AREA ESTIMATES - 1976 17 

3 - GENERAL CONCLUSIONS 20 

BIBLIOGRAPHY 21 


gfUgqgDIHQ £AGE BLANK NOT FILMK 


iii 



LIST OF FIGURES 


Page 


1 - The Jardinopolis Test Area 4 

2 - Visual interpretation of LANDSAT MSS imagery taken on May 10, 

1975 (Jardinopolis County) 13 

3 - Visual interpretation of LANDSAT MSS imagery taken on May 10, 

1975 (Jardinopolis County) 14 

4 - Visual interpretation of LANDSAT h6S imagery taken on May 10, 

1975 (Jardinopolis County) 15 


iv 



LIST OF TABLES 


Page 


2 

1 - Bare Soil area estimation (Km ) at different scales by four 

interpreters 6 

2 - Estimated arboreal vegetation area (Km ) at di fferent scales 

by four interpreters 7 

3 - Analysis of variance of bare soil estimated area 8 

4 - Analysis of variance of estimated arboreal vegetation area. 9 

5 - Landsat fteS Spectral responses of bare soil, arboreal vege- 

tation, and pasture/sugar cane for the analysed frames .... 12 


2 

6 - Preliminary results of areal crop estimates (Km ) in the Mu 


nicTpio of Jardinopolis by visual and automatic interpre- 
tation 16 

2 

7 - Estimates crop percentage and area (Km ) based on a sample 

of 63* CIR Images of Jardinopolis (1976) 18 

8 - Comparison of crop areal estimates using aerial photographs 

in both aircraft experiments (1975 and 1976) 19 


v 



ABSTRACT 


This paper' reports the results of several remote sensing 
techniques used for applied crop area estimations in an agricultural 
region of Sao Paulo State . Principal data sources included agricultural 
statistics from the Secretariat of Agriculture > CIR imagery acquired 
from t o INPE flight missions and multidate LANDSAT CCT's which were 
processed on the " Interactive Multispectral Imagery Analyser System" - 
INPE 's Image-100 . Efforts to produce area estimates concentrated on 
such ground features as "bare soil" (soil prepared for planting) and 
sugar cane. A total of 13 classes were sampled using CIR imagery and 
four classes - cotton , sugar cane , soybean and pasturelund - were 
interpreted from LANDSAT CCT's through the Image-100 processor. Flight 
experiments conducted in 1976 proved that the use of CIB imagery at a 
scale of 1:10 s 000 instead of 1:2.0,000 (1975), did not improve on the 
interpreters ' crop identification. An interpretation test to determine 
an optimum scale for the Image-100 classification of "bare soil " and 
"arboreal vegetation " showed no significant difference between scale 
changes but a scale of 1:150,000 was selected for training after 
considering field size within the study area. Areal estimates were made 
using a dot grid method and the Image-100 system. LANDSAT CCT's were 
partitioned into quadrangles and "single cell" and "clustering syntesis" 
algorithms were used for classification. "Bare soil" comprised 34.0% of 
the study area. A total of 156 Image- 100 hours were spent in the study. 



1 - INTRODUCTION AND APPROACH 

Brazil is essentially an agricultural country. It is the 
world’s largest coffee producer, second largest corn producer and a 
leading producer of sugar, cotton, rice, soybean and peanuts. The 
dynamics of cropping and food demand are subject to many variables which 
points to the necessity for timely crop production estimates. 

Considering that the current methods of crop prediction 
in Brazil are subjective and often outdated, the need for new techniques 
is essential in order to obtain greater accuracy and more timely 
information on a regular basis. The conventional system of crop forecasting 
which has been used for many years by the Secretariat of Agriculture, 
is based on the sampling of a number of agricultural properties where 
questionaires are administered 5 times during the growing season. It has 
been observed that the main source of error is in the area estimates 
given by farmers, while estimates of crop productivity are more accurate. 
For Brazil, LANDSAT-f-KS data processing offers new perspectives to meet 
our goals with respect to crop surveys over large regions. 

The Brazilian Institute of Space Research (INPE), in 
agreement with the Secretariat of Agriculture of Sao Paulo State, has been 
carrying out an Agricultural Statistics Project (EAGRI) since 1975. 

The objective of this project is to verify the feasibil ity of remote sensing 
techniques for crop forecasting, primarily using LANDSAT data and the 
"Interactive Mul tispectral Imagery Analyser System" (IMAGE-100). 

According to the project objectives INPE has had the 
responsibility of developing a metnodology using remote sensing techniques 
for crop area measurement. The Secretariat of Agricul ture has contributed 
ground data and productivity estimates. 

There. are many complexities inherent to a project of this 
type. Problems encountered during the development of a remote sensing 
methodology included: 



1) the growing season of annual crops is from October to May, in 
coincidence with the rainy season which restricts the aquisition 
of cloud-free imagery during this period of interest; 

2) major crops are planted in relatively small areas. This causes 
some difficulties in either spectral discrimination or in area 
estimates due to the low resolution of LANDSAT imagery. 

Considering these problems, and the current potentialities 
of remote sensing techniques, efforts to produce area estimates ‘ 
concentrated on the following ground features: 

1) BARE SOIL (area estimates of soil being prepared for planting). 
This information is important as it can be associated with the 
Secretariat of Agricul tures ' s publication of "The Report of 
Planting Intentions of Farmers". Both, area estimates and the 
Secretariat' s publications, could provide a preliminary crop 
prediction at the beginning of the growing season. The soil is 
prepared for planting from June to October. This is the best 
period of the year to acquire cloud-free coverage. Also, during 
these months bare soil is well defined on LANDSAT imagery due 
to the sharp contrast with neighboring vegetation. 

2) SUGAR CANE - area estimates of sugar cane can also be carried 
out during the dry season. Another advantage is that sugar cane 
is cultivated over a large and continuous area, which helps in 
crop identification. 


3) OTHER CROPS - the feasibility of using remote sens.ing data for 
coffee, citrus and pastureland surveys is also being studied. 


2 - INTERPRETATION METHODOLOGY AND RESULTS 


The study area (Figure 1) is the Regional Administrative 
‘vision of Ribeirao Prato (DIRA-RP). This region includes our test area, 
the MunicTpio of Jardinopolis, which contains the major representative 
crops of Sao Paulo State. 

2.1 - VISUAL INTERPRETATION 

Two CIR aerial surveys were conducted over Jardinopolis. 

The first was in March, 1975 (scale 1:20.000) and the second was in 
February, 1976 (scale 1:10.000). The objectives of both experiments were 
to test the feasibility of using aircraft data for area estimates and 
to give support to LANDSAT data interpretation. 

The EAGRI Project frequently uses areal measurements taken 
from CIR imagery to support LANDSAT data analysis. Two different methods 
of area measurements were compared: the dot grid method (using 1cm 2 and 
2cm 2 grids) and the weighing method. 

Five aerial photographs from 1975 were analyzed. The 
maximum and minimum percentages of ground cover for interpreted crops 
ranged from 50% to ]%, From the results of the statistical analysis, no 
differences were observed between the two methods. However, the dot grid 
method is more time-saving (using computer analysis) than the "weighing" 
method and thus the 2 cm 2 grid was selected for the crop area estimates 
in our study. 


The total area of Jardinopolis (521 km 2 ) was mapped by 
visual photointerpretation and supported by ground observations. A thematic 
map was constructed to show the boundaries of each crop present in the 
test area. Area estimates were made for each crop type. Both the thematic 
map and the area estimates were the first basic data compiled of the test 
site from IMAGE-100 interpretation. 




2.2 - AUTOMATIC INTERPRETATION 


An optimum scale determination for the IMAGE-100 operation 
was conducted to produce an accurate and time-saving classification of 
•'bare soil" and "arboreal vegetation". An area of approximately 
1,400 km 2 was selected and analysed at four different scales: 1:125,000, 
1:150,000, 1:200,000 and 1:300,000. A computer-aided interpretation 
(IMAGE-100) was administered by four trained technicians who were asked 
to classify bare soil and arboreal vegetation in four scales for one - 
third of the study area. A two way analysis of variance was performed. 

The interpreter's area estimations were run against the four different 
scales (tables 1 and 2). 

The results show that there is a significant difference 
among interpreter's scores for bare soil and arboreal vegetation at 
p = 0.10 and p = 0.05, respectively (tables 3 and 4). 

The results of tables 1 and 2 show that the individual 
interpretations of arboreal vegetation were more consistent than for bare 
soil at different scales. One can conclude that the interpreters should 
have a better a priori knowledge of the area. In the case of bare soil, 
the lack of field information in areas with different soil types and 
the difficulty of perceiving pixel density explained the differences 
among the interpreters ' classification performances. Arboreal vegetation 
is a class which can be easily identified without field information. 

The variation of pixel density is the prominent factor in explaining 
the differences among the interpreters' scores. This suggests that, 
besides using training samples, a proficient knowledge of the test area 
is required to quantitatively evaluate the capabilities of the IMAGE-100 
System. The results showed no statistically significant differences 
among the scores between scale changes. 

Considering the general field size of the region and the 
selection of representative training samples, the scale of 1:150,000 proved 
to be the most suitable for computer-aided interpretation and analysis. 



TABLE 2 



8,29 

















0.05 = 3 72 


2,3 - AREAL EVALUATION OF BARE SOIL, ARBOREAL VEGETATION AND PASTURE/ 
SUGAR CANE USING LANDSAT AUTOMATIC INTERPRETATION 


After an optimum scale was determined (1:150,000), a bare 
soil estimation was conducted for areas of the DIRA-RP region using the 
I MAGE -100 System. The following corresponding computer compatible tapes 
(CCT) were analysed: 


ORBIT/POINT 

DATES 

178/26 

09/13/75 

178/27 

09/13/75, 05/10/75 

192/26 

09/14/75 

192/27 

09/14/75 


Support data included four LANDSAT MSS frames (channels 5 
and 7) at scales of 1:1,000,000 and 1:250,000, CIR imagery at scales of 
1:10,000 and 1:20,000 and ground information. 

Each frame of the whole scene was partitioned into 
quadrangles equival ent' to l/-|g of the original LANDSAT frame. This provided 
a scale of approximately 1:1 50,000 on the Image-100 display. 

Initially, several algorithms of the Image-100 System were 
used. Then, spectral responses of the targets (i.e. bare soil, arboreal 
vegetation and pasture/sugar cane in the four MSS channels were acquired.) 
This preliminary analysis led to the selection of "single cell aquisition" 
and "clustering synthesis" algorithms for classification purposes. 




Quadrangles were classified according to these algorithms. The spectral 
resolution used was 128 bits because a 256 resolution did not improve 
the classification and a 64 resolution generated too much overlapping of 
themes. 

The results show that 12,678 km 2 were classified as bare 
soil for September 13, 1975, This corresponds to 34% of the D1RA-RP 
region which has an area of 36,801 km 2 , as assessed by the Image-100. 
Besides bare soil evaluation, pasture/sugar cane (as a single class) 
and arboreal vegetation were also classified. These three classes, 
together, gave an area of 30,637 km 2 which corresponds to 83% of the 
DIRA-RP area. A pixel by pixel analysis will be carried out in the 
future, to determine classification accuracy by comparing the automatic 
classification output with aerial photographs. A total of 156 Image-100 
hours were spent in this study. Only 75 hours were effectively spent 
for bare soil classification, and the rest was spent for radiometric 
correction and recording. 

Table 5 shows the spectral responses in the four channels 
of the four LANDSAT frames analysed for bare soil, arboreal vegetation 
and pasture/sugar cane. These results are the average of the responses 
(upper and lower bounds) of each class, in all quadrangles of each 
frame mentioned previously. 

2.4 - CROP IDENTIFICATIONS AND AREA ESTIMATES - 1975 


The preliminary results of the 1975 crop inventory work 
for the Jardinopol is test area (Fig. 1) are discussed below: 

Figures 2, 3 and 4 show area delimitations visually 
interpreted from LANDSAT channels 5 and 7 (May 10, 1975) at a scale of 
1:250,000. This particular scene was concomitant with the harvest of the 
major crops. The results of areal crop estimates, based on visual 
interpretation of UR imagery and automatic interpretation of LANDSAT 
CCT's, are presented in Table 6. 


LANDS AT. MSS SPECTRAL RESPONSES OF BARE SOIL , ARBOREAL VEGETATION 










































Visual interpretation of LANDSAT MSS imagery taken on May 10, 1975 ( Jardinopol is County) 



TABLE 6 


PRELIMINARY RESULTS OF AREAL CROP ESTIMATES (KM 2 ) IN THE MUNlClPIO O F 
JARDIN0POLIS BY VISUAL AND AUTOMATIC INTERPRETATION 


CLASSES 

" ■ •' ’ ,u - ■ " T 

AREAL ESTIMATE (KM 2 ) IN 1975 

PERCENT 

CHANGE 

B - Ax 100 
A 

VISUAL INTERPRETATION 
USING CIR FILM V A) 

AUTOMATIC INTERPRETATION 
OF LANDSAT DATA (B)* 

COTTON 

42.31 

33.40 

- 21 .06 

SUGAR CANE 

25.09 

21.40 

- 14.71 

SOY BEAN 

66.06 

70.67 

+ 6.97 

PASTURELANC 

182.86 

166.13 

- 9.15 
— ■ . - 


* LANDSAT CCT May 10, 1975. 

Comparing the percent change between the four classes, 
soybean was overestimated on the IMAGE-100 compared to the conventional 
interpretation method (+ 6.97%). A moderate information loss was scored 
for pastureland {- 9.15%) and increasing losses of interpretabil i ty were 
shown for sugar cane (- 14.71%) and cotton (-21.06%), which ranked last 
among the four classes. It is apparent that automatic interpretation 
more closely approximates visual interpretation with increasing size of 
the target area. Both the crop area delimitations and area estimations 
could be improved if CIR imagery and sequential LANDSAT data were 
acquired during the growing season. 




2,5 - CROP IDENTIFICATIONS AND AREA ESTIMATES - 1976 


In the 1976 flight mission over the Jardinopol is test area, a 
scale change to 1:10,000 was made with the objective of determining if 
greater accuracy of crop identification could be achieved. The aerial 
survey was conducted during the first fortnight of February, midway 
through the growing season, for the major crops. The visual interpretation 
of the imagery was supported by ground truth data, which combined to give 
100% identification accuracies for: cotton, coffee, sugar cane, pastureland, 
soybean, orchard and cerrado. High values were obtained for corn (98%) 
rice (94%), peanuts (92%), and sorghum (82%). The use of CIR imagery at 
a scale of 1:10,000 (1976), instead of 1:20,000 (1975), did not improve on 
the interpreters' crop identification accuracy. It is concluded that 
with the high identification accuracies achieved at a scale of 1:20,000, 
smaller scales could be used to reduce operational and interpretation 
costs. 


Table 7 presents estimated crop percentages and equivalent 

2 

areas (km ) from a sample of 63 CIR images. Errors of estimations were 
calculated for 13 classes. The mean percent error of estimation was 
21.69%. A total of 100 images cover the Jardinopolis area. A calculation 
was made to determine the number of photos needed for crop area estimations 
with a 10% error. After the termination of this phase of the project it 
was concluded that the CIR imagery at a scale of 1:10,000 was inadequate 
for crop inventory work. First, there was no improvement in crop 
identification accuracy compared to a scale of 1:20,000. Second, the 
photographic sampling units at this scale produced a high variance 
in the crop estimations. 

Table 8 shows the results derived in both the 1975 and 
1976 aerial CIR surveys. 

The comparison of these two tables made possible the 
verification of crop changes from one year to another. For example, it 
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TABLE 8 


COMPARISON OF CROP AREAL ESTIMATES USING AERIAL PHOTOGRAPHS IN 


BOTH AIRCRAFT EXPERIMENTS (1975 AND 1976) 


CLASSES 

2 

ESTIMATED AREAS (KM*) 

PERCENT DIFFERENCE 

1975 

1976 

BASED ON 1975 DATA 

Cotton 

42.32 

39.65 

- 18.12 

Peanut 

4.52 

5.00 

+ 10.62 

Rice 

4.76 

11.67 

+145.17 

Coffee 

7.94 

5.84 

- 35.96 

Sugar Cane 

• 25.09 

52.37 

+108.73 

Corn 

81.38 

75.66 

- 7.03 

Soy bean 

66.06 

40.59 

- 38.56 

Orchard 

31.55 

21.88 

- 30.65 

Pastureland 

182.86 

178.83 

- 2.20 

Cerrado 

35.19 

27.41 

- 22.11 

Artificial Forest 

2.53 

2.61 

+ 1.16 

Sorghum 

- 

10.06 

- 

Other 

36.82 

54,50 

+ 48.02 

TOTAL 

521.07 

521 ,07 

, 

- 










was observed that areas planted with sugar cane increased more than 
100%. This induced a decrease in areas planted with other crops. 

3 - GENERAL CONCLUSIONS 


The following conclusions were drawn from this stage of 
the EAGRI research and development project: 

1} To meet the objectives of crop forecasting using remote sensing 
techniques, it is essential to observe the spectral changes of 
the target crops during their growing season using sequential 
imagery. 

2) The two flight experiments (1975 and 1976) showed that an 
effective scale change from 1:20,000 to 1:10,000 did not improve 
the accuracy of crop identification. Consequently, the increase 
in the quantity of imagery at 1:10,000 only increases time and 
operational costs. 

3) Quantitative techniques and digital processing of LANDSAT data 
in the IMAGE-100 System permitted rapid automatic interpretation 
and areal estimations. These results were encouraging and will 
lead to the increasing use of the IMAGE-100 System for the 
analysis of large agricultural regions. 


